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Abstract 

Within the context of the functional renormahzation group flow of gravity, we suggest that a 
generic f{R) ansatz (i.e. not truncated to any specific form, polynomial or not) for the effective 
action plays a role analogous to the local potential approximation (LPA) in scalar field theory. 
In the same spirit of the LPA, we derive and study an ordinary differential equation for f{R) 
to be satisfied by a fixed point of the renormahzation group flow. As a first step in trying to 
assess the existence of global solutions (i.e. true fixed point) for such equation, we investigate 
here the properties of its solutions by a comparison of various series expansions and numerical 
integrations. In particular, we study the analyticity conditions required because of the presence 
of fixed singularities in the equation, and we develop an expansion of the solutions for large R up 
to order N = 29. Studying the convergence of the fixed points of the truncated solutions with 
respect to A^, we find a characteristic pattern for the location of the fixed points in the complex 
plane, with one point stemming out for its stability. Finally, we establish that if a non-Gaussian 
fixed point exists within the full f{R) approximation, it corresponds to an R^ theory. 
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1 Introduction 

The application of functional renormalization group techniques to gravity has generated many 
appealing results in support of the asymptotic safety scenario [H [21 [3l HI [5] . The main tool used 
in such investigations is the so-called Functional Renormalization Group Equation (FRGE) for 
the effective average action [6], which read^ 



4-rfc[$] = -STr 
dt ^ ^ 2 



I ^ + TZk -TZk 



(1.1) 



The FRGE is an exact equation whose solutions determine a flow of effective actions in the theory 
space of all possible functionals F^, interpolating between a bare action at some initial UV scale 



^Here $ denotes the collection of all the fields in the theory under consideration, and STr a functional supertrace 
over their spinorial indices and spacetime coordinates (collectively denoted by the A, B indices). The running scale 
hit = Infc, and TZk is a cutoff function implementing the Wilsonian momentum-shell integration. For further details 
we refer to the many general reviews [71 El El IIOI llll I12| . 
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k = A and the full effective action at A; = 0. The primary goal of the asymptotic safety program 
is to show that there exists a non-trivial (i.e. non-free) gravitational action T* which is a fixed 
poin1§| for such a flow, with finitely many relevant directions. 

Due to the complexity of the question, the main line of progress in this topic has been based 
on a sort of "mathematical experiments": an equation, which is in principle exact, is solved 
by truncating its infinite-dimensional functional domain to a finite-dimensional subspace; such 
operation is repeated for various truncations, and stability and convergence of the results are 
tested. In practice, exploring larger and larger truncations is a very tough and tedious job, and 
the question is how confident can we be about our conclusions, based on the results we have 
obtained so far. 

Such a question is of course not specific to gravity, but common to any application of exact 
renormalization group equations in which the approximation scheme does not rely on a small 
parameter expansion, but rather on an apparently arbitrary truncation of the space of action 
functionals. In simpler settings than gravity, as for scalar field theory, the use of the FRGE 
has reached a sufficient level of confidence that allows us to make very solid statements about 
the phase structure of the theory and even quantitative predictions that can compete with other 
methods (see for example [13] and references therein) . One essential element for such achievement 
is the use of the RG equations as differential equations for unspecified functionals, rather than 
for the couplings of a usual field expansion. For instance, rather than truncating the scalar 
potential Vk{(j)) to a polynomial of order A^, looking at the RG flow of the coefficients of the 
various monomials and studying what happens as N is increased, one can instead view the RG 
flow equation as a partial differential equation for the unknown function Vk{(f>)- The lowest order 
of approximation in such a scheme is called Local Potential Approximation (LPA) and consists in 
retaining only the equation for 14 (</>). At next order one includes the equation for a "wave function 
renormalization functional" Zk{(j)), and so on. Explicitly, one writes a derivative expansion of the 
type 



which plugged into (jl.ip leads to partial differential equations for the unknown functions 14 ((/)), 
Zfc(0), etc. A comparative study between such scheme and the polynomial truncations was 
carried out by Morris in [13], highlighting the greater reliability of the former, which leads to 
more accurate results, and to a better understanding of the convergence properties of truncations. 
Of course polynomial truncations are still used because they are easier to handle, and they can 
provide accurate results in many cases (e.g. [15]), however the validity of such truncations is 
better assessed by a comparison to the derivative expansion. 

^That is, dtr* — 0. Here dt is a partial derivative, whereas in {O} d/dt is a total derivative, the difTerence 
being that the latter acts also on the scaling dimension of the fields and of the Lagrangian itself. This point will 
be made more explicit in (|2.5p . 




J d^x[Vk{cl>) + Zki<p)d''<pdf,ci) 



(1.2) 
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In the case of pure gravity most of the truncations studied so far take the forn: 

TV 



i=0 

where R is the Ricci scalar, and the truncation order has been increased from the original = 1 
[m dZl [m HI], and iV = 2[20l|2l],toAf = 6[22l|23],Ar = 8[2l]andAf = 10 [25]. In t26j an 
RfjiupcjR^'^'"' term was added also to the N = 2 truncation, and in [23] some non-analytic terms 
like R~^ and Ini? have been included]^ 

Organizing the truncations of the gravitational action in powers of curvature is a natural 
starting point, and one familiar in effective field theory [33]. Such truncations do not correspond 
of course to polynomial truncations in the fundamental field, which is the metric, and which 
appears non-polynomially in the curvature, thus differing in this respect from the polynomial 
truncations in the scalar field theory case. On the other hand, just as in that case, polynomial 
truncations lead to ordinary differential equations for the flow of a finite number of couplings, and 
to algebraic equations for the fixed points. Furthermore, not all higher-order invariants in the 
curvature contribute with higher-order derivatives to the two-point function, thus the expansion 
in powers of the curvature is in this sense also not a derivative expansion. 

In the present work we advocate the point of view that truncations of the type p.3p are akin 
to the polynomial truncations in scalar field theory, and we propose an analogy between a generic 
f{R) approximation and the LPA. Of course, due to diffeomorphism invariance, no potential 
can be written just for the metric fluctuations, and instead the simplest Lagrangian that can 
be written without restricting to any specific function is that of an f{R) theory. That is, we 
retain a generic dependence of the action on the scalar curvature R, but discard any derivatives 
of R as well as any more complicated tensorial structures like R^yp^R^'^^'^ ■, etc. It is well known 
that such a theory has only one additional (scalar) degree of freedom with respect to general 
relativity [341 135j. or in other words it only adds fourth-order derivatives to the trace component 
of the metric fluctuations, which in general relativity is not a dynamical degree of freedom. 
Fourth-order derivatives for the transverse-traceless sector of the metric would be contained for 
example in terms like f2{RpupaR^'^^'^), while higher derivatives for both sectors would start for 
example with terms hke Zi{R)RV^R or Z2{R)R^'^V'^Rpy. We will not attempt here to define 
a full derivative expansion to higher orders, however it is clear that the f{R) functional is the 
functional ansatz with the minimal number of derivatives (of course among those containing an 
arbitray function, Einstein-Hilbert or Gauss-Bonnet being special cases), and hence it can serve 
as the leading term in such an approximation scheme. 

There is also another analogy that can be drawn between the LPA and the f{R) ansatz. In the 
evaluation of one first expands an action functional like p.2|) around a background field i^, 
writing <j) = (p+ip and computing the second variation with respect to 99, then one plugs the result 



^Details of the contruction of truncations in the gravitational case will be reviewed in Sec. [2] 
^Matter has also been included [271128) . as well as corrections to the ghost sector |29II30| and bimetric truncations 
[311 132| . Here we will concentrate on single-metric pure gravity sector of the theory. 
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into (|l.ip and projects it on the background, setting ip = 0. In light of that, we can restate the LPA 
as being the approximation in which we take a constant background (p: once we project the FRGE 
on such a background, clearly the only functional we can discern is a potential. Analogously, in 
the gravitational case once we choose a maximally symmetric background (typically of spherical 
topology), as it is usually done in order to keep the calculations manageable, the only functional 
we can discern is an f{R). Indeed for a spherical background the Weyl tensor is identically zero, 
and the Riemann and Ricci tensors are both proportional to the Ricci scalar R, which is constant. 
Note that given the approximation consisting in the choice of a spherical background, no further 
approximations are needed. In particular no truncation is needed, as we can write down a generic 
action which contains all possible terms that might ever be generated by the functional traces of 
the FRGE on a spherical background: this is the f{R) action. 

In the present work we will study the gravitational FRGE in the f{R) approximation, in the 
spirit of the LPA. We will derive the ordinary differential equation to be satisfied by a fixed- 
point f{R) function, and we will study general properties of its solutions. As pointed out in 
|36j for the LPA and emphasized in [TJ], most of the solutions to the fixed-point equation end 
at a singularity, while a putative fixed point should be represented by a non-singular solution^ 
For our equation the identification of global solutions turns out to be a very challenging task, 
and we report here the present status of our understanding, postponing to future work a more 
comprehensive numerical study of the solutions. Here we will examine in some detail the different 
type of singularities that plague our equation, identify two analyticity conditions to be satisfied 
by the solutions, and develop a new series expansion, for large R, which appears to be much more 
manageable than the usual expansion at i? = 0. The main practical outcomes of our work are an 
identification of candidate fixed points from the expansion of the solutions at large R, and the 
observation that if a non-trivial fixed point exists, it necessarily corresponds to an R^ theory. 

In Sec. [2] we present more precisely the f{R) ansatz, and discuss the main results of the paper. 
The reader who is not interested in the technical details of the computation can in principle skip 
directly to the conclusions after reading Sec. El All the technical details will be presented in 
the remaining sections: in Sec. [3] we provide the functional variations of the action functional, 
together with the ghost and auxiliary sectors; in Sec. H] we present the cutoff scheme employed 
in our computation; in Sec. [5] we describe the method used in evaluating the functional traces. 
Finally in Sec. Owe collect the results into the final form of our FRGE, and in Sec. [7]we detail the 
analysis of the fixed-point equation. A summary of results and future prospects ends the paper 
in Sec. [HI 

^Besides being a reasonable physical requirement, it was also proved by Felder [57], within the LPA, that any 
fixed point which is the limit for t — > oo of an effective potential with non-singular initial condition at t = must 
be a global solution of the fixed-point equation. 



5 



2 The f{R) approximation: setup and outline of results 

The flow equation (jl.ip is adapted to the gravitational case along the lines of the standard fleld- 
theoretic quantization, as used also in one- and two-loop calculations |38| I39|. Difi'erences only 
reside in the cutoff choice, and in the approximation being used. 

In the case of pure gravity, the fields comprise the metric, the ghosts and occasionally some 
auxiliary fields implementing the functional Jacobians originated by field redefinitions. We use 
the background field method to obtain a gauge-invariant average effective action, and we define 
the decomposition of the metric by 

l^,u = 3^lu + h^y , (2.1) 

with g^jy denoting the background and h^^ the fluctuations. As background metric we will take 
a d-dimensional sphere, as explained in the introductionlfl 

Following |16] . it is useful to cast a general truncation of the effective average action in the 
following form: 

Tfcit, $] = fk[g\ + f g] + Tgf [/i, g] + rgh[/i, g, ghosts] + 5auxb, aux.fields] . (2.2) 

In this decomposition Tk{g] depends only on the total metric, and it is the proper gravitational 
action. Tgf and Fgh denote the gauge-fixing and ghost-terms respectively, for which we will take 
the classical functionals but eventually allowing a running of the gravitational couplings, while 
5aux is a coupling-independent quadratic action encoding the Jacobians. Tk[h,g\ vanishes for 
/i = 0, and it encodes the deviations from standard Ward identities due to the use of a cutoff [16] . 
The role of such term has been investigated via bimetric truncations in [31^ [32] . In the present 
work we will make use of the common approximation = 0; in such case, it suffices to study 
the FRGE at h^^ = 0. 

Our ansatz for the gravitational action is 

fk = Zk j d'^x^hiR) , (2.3) 

where = (IGvrGfc)"^ and is the (running) Newton's constant. The running RG scale is fc, 
and we will introduce the dimensionless Ricci scalar R = Rjl? ^ and Lagrangian 

whose shape we will try to fix by use of the FRGE. 

When plugging our ansatz (j2.2ll2.3p into p.ip we will obtain a partial differential equation for 
fk{R)- We can trivially write down the left-hand side of the FRGE: 

d 



dt 



I h^u —ghosts— 



j d'^x^ [dtUR) + dfk{R) - 2Rfl{R)] . (2.5) 



®For what follows, all one needs to know about the sphere is that the Ricci scalar is constant, i?^^ = ig^^^i?, 
Rfivpa = d{d-i) (ffMpffi^g- ^ gficTgvp)R, and that the radius and the volume are given by = d{d — 1)/R and 
V = {4:Tvp^)'^^^r{d/2)/r(d), respectively. We will also need the eigenvalues of the Laplacian, reported in Table [1] 
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Note that dt here is a partial derivative acting only on the explicit /s-dependence in the function 
fk{R), while primes denote differentiation with respect to R. The right-hand side of the equation 
is also a function of f{R) and its derivative with respect to both k and R. Hence the FRGE is a 
partial differential equation (PDE) for the function f{R)- 

Now remember that a step of the RG flow consists of two intermediate steps [8] : integration 
over a momentum shell, which we take care of via the FRGE; and rescaling of all dimensionful 
quantities to restore the original UV cutoff. As we have translated everything into dimensionless 
quantities, the second step is trivially taken care of, and hence a fixed point of the RG flow has 
to satisfy 

dtfk{R) = 0. (2.6) 

As a result, the PDE obtained from the FRGE will reduce to an ordinary differential equation 
(ODE) at the fixed point. In contrast, in polynomial truncations one gets a system of ordinary 
differential equations for the flow, and a system of algebraic equations for the fixed point. As 
explained in the intro, the price we have to pay for going to a higher level of equations should be 
compensated by the gain in confidence in the results thus obtained. 

We are now going to highlight the main results of our analysis, and anticipate some conclusions. 
All the details will be given in the following sections. 

The study of the FRGE in the described approximation presents two challenges: the actual 
evaluation of the r.h.s. of the equation (|l.ip . and the numerical study of the PDE obtained in 
this way. The first part will be described in sections O HI [5] and [6l A similar calculation has 
been done before \23\ 124]. but we will adopt here a different implementation of the ghost sector, a 
different cutoff scheme and a different evaluation technique for the traces, mostly following |40| . 
The main result of this part of our work is hence is the derivation of the fixed-point equation in 
d = 4 dimensions, which takes the fornj^ 

f"iR) = M{f,f',f";R) ^ ^ ^2.7) 

i2(i?4 _ 54i?2 _ 54) hp> _ 2)f'{R) - 2f{R)\ 

where AA(/, /', /"; R) is a polynomial in f{R) and its first two derivatives, with coefRcients poly- 
nomial in R. Its precise expression will be given later in (|6.7p . 

For the second part (Sec. [7]), we will restrict to the question of existence of fixed points, thus 
reducing the task to the study of the ordinary differential equation (|2.7p . We will address in 
particular the following topics: 

1. fixed singularities, initial conditions and analyticity conditions; 

2. movable singularities; 

3. large- i? expansion (boundary condition at infinity); 

4. fixed points from truncations of the large- i? expansion; 

5. generic form of a fixed-point effective action. 

^We denote the function at the fixed point just by omitting the subscript k from /. 
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As the solution of a third-order ODE is specified by three initial conditions, it would seem 
at first that one would obtain from (|2.7p a continuous of solutions parametrized by the initial 
conditions. In particular, denoting the initial conditions a^ 

/(0)=<.o.5i.. /'(0) = <.,.-^, n„) = <.,.^, P.B) 

it would seem that the cosmological constant A, the Newton's constant G and the higher-derivative 
coupling c* are completely free at the fixed point. However, an analysis similar to [H] should 
lead to the conclusion that only a finite subset of solutions satisfy basic regularity requirements, 
like not having singularities at finite R. Indeed we find that for generic initial conditions the 
solution develops a logarithmic singularity (see (|7.7p ) at a finite value of R. We distinguish two 
type of singularities, movable singularities and fixed singularities. Singularities of the first kind 
appear because of the non-linear nature of the ODE, and as it is suggested by the name, they 
occur at a location that varies with the initial conditions. On the contrary, fixed singularities 
occur at those values of R where the equation is explicitly singular. In the specific case, the latter 
correspond to the zeros of the denominator in (|2.7|) . which has three real zeros, at i? = and at 



R = ±^3(9 + v87) = R±. The presence of such fixed singularities will play an important role 
in our analysis. 

A brief explanation is due also for the fact that the equation (j2.7p is third-order, given that in 
(jl.ip only the second functional derivative of F appears. The reason is that, following a standard 
procedure [231 111], we adapt the cutoff TZ^ to the Hessian of the action (see Sec. H]), meaning 
that the cutoff contains the second derivative fl!{R). When the total derivative d/dt acts on the 
cutoff in (jl.ip . we obtain a term — 2ii/^"(i?), among others. This explains at the same time why 
the equation is third-order, and why it has a singularity at -R = 0: that is a point at which the 
coefficient of f"'{R) vanishes. The appearance of such a fixed singularity is hence very generic, and 
it actually acts in favor of a qualitative scheme-independence: it is well known that a quantitative 
analysis of the FRGE depends on the cutoff choice, but one would expect that there would be 
even a qualitative difference between a third order (in the present scheme) and a second order 
(in an hypothetical scheme with no f'^{R) dependence in the cuttoff) equation, i.e. a different 
number of initial conditions could lead to a continuos of fixed points or no fixed points in one 
or the other case. However, in presence of a singularity for the third-order equation the number 
of independent initial conditions at i? = is reduced to two by the requirement of regularity of 
the solutiorI§: upon substitution of a Taylor expansion f{R) = X^„>o c-nR" we find that to lowest 



^Obviously the choice of parametrization of the initial conditions is driven by the identification of the first two 
terms with the corresponding ones in the Einstein-Hilbert action, with A and G the dimensionless cosmological 
and Newton's constants. The star denotes the fixed-point values of such parameters. 

®A similar situation occurs in scalar field theories [41] if as fundamental field one uses y = (f^ instead of (f). 
In that case the regularity condition replaces the even potential condition V {(f) = 0) = 0. We verified that the 
same could be done here, defining U{x) = f{x^), and the regularity condition being traded for the initial condition 
!7'(0) = 0. As nothing is gained in such scheme, while the gravitational interpretation gets a bit obscured, we stick 
in the following to the variable 7?. 



8 



order in R the equation (j2.7p imposes a condition, 



S_(ao, ai, 02) = 24ao + 17aiao + 216a2ao 

+ 3847r^ (ao + ai) (2ao + 3ai + 1802) ao - 27a? - 144aia2 = 



(2.9) 



which for example can be easily solved for 02. 

The reason for the singularity at Rj- is quite technical and will be explained in some detail 
in Sec. 17.21 Roughly speaking it has to do with the constant scalar mode. We found its nature 
to be also quite generic, and we could not devise any reasonable scheme to make it disappear. 
On the contrary, by a clever choice of cutoff scheme we were able to eliminate a number of other 
fixed singularities which appeared in previous versions of the equation |23^ I24j. Similarly to the 
singularity at the origin, we identify an analyticity condition at R±, equation (j7.4p . relating the 
first three coefficients in a Taylor expansion of the solution at the singular point. 

Unlike for finite R, solutions are allowed to diverge at infinity. We find that a solution for R — t- 
00 must behave like f{R) ~ AR^, with sub-leading corrections (detailed in Sec ]7.4p dependending 
on the single free parameter A. Such single parameter dependence provides an easier setting for 
numerical investigations (as compared to i? = where we have two free parameters), which can 
been performed by shooting backward from given initial conditions at large R. 

A true fixed point should correspond to a global solution, meaning a solution satisfying the 
analyticity conditions at the fixed singularities, presenting no other singularities at finite R, and 
matching the large R expansion at infinity. Due to the complexity of the equation, the quest for 
such a global solution turns out to be quite challenging, however with the help of preliminary 
results from numerical integrations we will show that these seemingly too many conditions can 
in principle be satisfied. 

Postponing a comprehensive numerical investigation to future work, here we will exploit fur- 
ther the new insights that can be derived from the large R expansion. In particular we find that 
treating such expansion as a standard truncation, we can identify an interesting structure for 
the fixed-point solutions (Fig. [3]), and we can single out one point with surprising convergence 
properties. Unfortunately at such point the number of relevant directions seems to be increasing 
with the order of the truncation. 

Finally, from the same large- i? asymptotic expansion, we conclude that if a global solution 
to the fixed-point equation exists, it must necessarily correspond to a fixed-point theory with 
effective action F* = A* f d'^x^/gR^ , for some finite A* . Such a result is also compatible with 
recent works [42^ H3| , in which the identification /c^ ~ i? is assumed for truncated versions of 
the average effective action. Of course an R^ action could also be expected on simple dimensional 
grounds, as it is the only f{R) action possessing scale invariance, however it had never been 
verified before that this is indeed the case, as it is not obvious to see such scaling from truncated 
expansions, and it could in principle be not the case if the Ricci scalar had to acquire an anomalous 
dimension at the fixed point. Consistently with the LPA analogy we find here that the scaling is 
not anomalous. 

The following sections will provide all the details for the interested reader. 
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3 Hessian and gauge-fixing 



In order to evaluate the r.h.s. of the FRGE, the first step is to compute the second variation of 
the ansatz at hf^i, = 0, i.e. the Hessian on the background. 

Although functional variations of the action functional for f{R) theories have appeared before 
in the literature [23\ I24j . we report them here again for completeness, and because of a slight 
difference in the choice of variables and notation. 

We use the transverse-traceless decomposition of the metric fluctuations, given by 



h 



fJ,U 



r 1 - 



with the component fields satisfying 



V^C^ = 0, h = h + ^a, h = g^^h^"" 



(3.1) 



(3.2) 



and A = -V^. 

The Hessian for our ansatz comprises the following components: 



p(2) 



2 



R 



a 



d{d-l)' 



R 



R] + ifk- -Mfk 



2R „, 



^(2) 
hh 



d-2 



2{d-lf 
d-2 



/HA 



R 



d-2 



a 



2{d-l) 
d-2 



s(2) 
ha 



A 



d-l 
d-2 



+ id-l)fUA 



R 



d-l 



Rf'k - 2^ 



Zk- 



2(i2 



Rfk-^fk]A, 



R 
d-l 



+ a 



d 



d-2 



Rf'k - 



A 



R 



d-l 



A . 



We have already included in the Hessian the terms coming from the gauge-fixing action 

r^f[h,g] = ^ j d'^x^F^F'^, 



with 



V, 



d-l R\ 

(2) 



(3.3) 
(3.4) 



(3.5) 
(3.6) 

(3.7) 



(3.8) 



(3.9) 



In the gauge a — )■ 0, that we will use in the following, neither F^^ nor the other terms 

(2) (2) 

proportional to the equations of motion in Taa and T^^ contribute to the traces. 
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3.1 Ghosts and auxiliary fields 

For the ghost action we follow [40] and write 



(3.10) 



where the Cj and c are complex Grassmann fields, while BJ^ and b are real fields, and the index T 
denotes transverse vectors. As explained in [50], while being formally equivalent to the standard 
implementation, this version of the ghost sector has the merit of realizing an exact cancellation 
(generally on shell, but also off shell in the q = gauge) between ghosts and pure-gauge degrees 
of freedom, thus ensuring gauge-independence of the on-shell effective action. We have explicitly 
verified that all the qualitative features emerging in the following analysis are found also with a 
standard version of the ghost action. 

The final ingredient of our truncation is the action for the auxiliary fields, introduced to 
take into account the Jacobian arising in the TT decomposition (|3.ip . The Jacobian for the 
gravitational sector leads to the auxiliary action 



(3.11) 



where the xji, ™d x ^^re complex Grassmann fields, while (^J and C are real fields. The Jacobian 
for the transverse decomposition of the ghost action is given by 

5aux-gh = / d'^x^^A^, (3.12) 



with 6 a real scalar field. 



4 Cutoff scheme 

We will use here a variation of the scheme introduced in [30] , which we will call "on-shell" Type 
II cutoff (in the spirit of the nomenclature of |24j). 

Denoting with rfc(z) some fixed cutoff profile function, we recall that the scheme introduced 
in [30] amounts to choosing the cutoff TZk in such a way to implement in the on-shell part of 
Hessian the rule 

A-n(^).A + .V(^). (4.1, 

In the course of the present work we realized that certain unphysical singularities appearing in the 
functional traces (see [IQ], but also [23] and [21]) actually appear from this unfortunate choice. 
Consider for example the typical operators appearing in the scalar part, Aq = A — ^j^. When 
using (|4.ip in combination with the optimized cutoff (|5.2p . we obtain for the modes below k'^, 

Ao^k^-^, (4.2) 
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which is of course zero, and hence not invertible, at R = d—1. In order to avoid such singularities, 
we win adopt here the following set of replacement rules: 



Ao^A 
Ai = 
A + 



R 



A 



d-l 
R 

~d 
2R 



P, 



(o)/Ao 
U2 



d{d - 1) 



,(2)/A2 

U2 



A, + .V.(^), 



/A2 



(4.3) 
(4.4) 
(4.5) 



and the profile function rfc(z) will be chosen later. 

Note that as we will work in the gauge a = the difference between on-shell and off-shell 
scheme is not essential, and it is basically only used as a motivation for including in the T]^^ j, 

(2) 

and operators certain potential terms and not others 
In practice, we have the following cutoff functions 



^ a 



A? 



7^f 



= Zk 
Zk 



d-2 
W 

(d-l)2 



2{d-l) 
d-2 



'fk 



U2 



d^ 



a 



((^^^°^( J)) ^ d-l (^^°^(fc2 



K]+{d-i)fkk'n{^ 

R 



Ag(Ao + 



d-l' 



(4.6) 
(4.7) 
(4.8) 
(4.9) 



and so on for ghosts and auxiliary fields. 



5 Spectral sums 

We are going to evaluate the traces by a direct spectral sum rather than in a heat kernel expansion. 
By spectral sum we mean that a generic trace will be evaluated as 

Tr,T^(A) = ^ Dn,sW{Xn,s) , (5.1) 

n 

where {A„^s} is the spectrum of eigenvalues of the Laplace-type operators (j4.3ti4.5p on spin-s fields, 
with the relative multiplicities {Dn^s}- The spectra and multiplicities are reported in Table [TJ 
and are obtained by appropriately shifting the standard one (for example see [15]). 

We have to be careful not to include fictitious modes in the sum. Remembering our decompo- 
sition for the metric fiuctuations (|3.ip , we see that we should exclude two sets of modes that give 

^°We recall that the crucial aspect of the "on-shell type" cutofT introduced in [40j was to avoid introducing a 
gauge-breaking cutoff term, that is a cutoff term for the gauge-variant field components and a with no ghost 
counterpart. In the gauge a — the fields and a only survive in the gauge-fixing part of the Hessian, and the 
effect of the regulator on such terms is correctly taken into account by (|3.f0|) in combination with (|4.3l4.4p . 
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Spin s 



Eigenvalue A^^^ 



n{n+d—l)~d 
d(d-l) 



Multiplicity D^, . 



(n+rf-2)! {2n+d~l) 

n!(d-l)! 

(n+d-3)! n(n+d-l)(2n+rf-l) 
(rf~2)!(n+l)! 
(n+d~3)! {d+l)(d~2)(n+d)(n~l)(2n+d- 



i?; n = 0,l, 



n{n+d~l)~d 



R;n= 1,2... 



n{n+d—l) 
d{d-l) 



R: n = 2,3. 



2{d~iy.{n+Tjr 



Table 1: Eigenvalues of the Laplace- type operators ()4.3ti4.5p on the d-sphere and their multiplic- 
ities 



no contribution to /i^jjy. First, we should exclude the Killing vectors, satisfying V + Vjy^^j = 0. 
Second, we should leave out also the constant scalar modes a = constant. A similar set of modes 
should be excluded also from the ghosts and auxiliary fields, as these are all fields introduced 
hand-in-hand with ^ and a. The only fields for which we retain all the modes are hj^^ and h. Note 
that, differently from |24t I23|. we do not exclude the scalar modes corresponding to conformal 
Killing vectors = V^o", i.e. those scalar modes satisfying V^V,yf7 = \g^vS/'^(T. It is indeed 
clear that in our decomposition (j3.ip such modes do contribute to h^y. This can be seen also 
from the point of view of the ghosts: the ghost modes should be in one-to-one correspondence 
with the modes of the gauge parameter (e^,e), and from C^g^y = V^e^ + Vj^e^ -|- 2V^V^e it is 
obvious that there is no reason to exclude the scalar modes e corresponding to conformal Killing 
vectors. As a consequence the tensor and vector sums will start at n = 2, while all the scalars 
sums will begin at n = 1, except for the h mode starting at n = 0. 
We choose to work with Litim's optimized cutoff j46j 

rk{z) = (1 - z)e{l - z) , (5.2) 

for which 

dtik^rk{A/k^)) = 2k^e{k^ - A) . (5.3) 

Its great technical advantage is that with it all the functions appearing in the FRGE have a 
numerator proportional to the step function, and hence the spectral sums are cut off at Ng = 
max{n G N : \n,s < A:^}. At the same time, for all An,s < A:^, we have Pki^n^s/k"^) = k'^- 
We write the FRGE as 

N2{R) Ni{B) 
dtrk= W2{Xn,2/k^R)+ Wi{Xn,l/k\R) 

n=2 n=2 

NoiR) mm (5.4) 

+ E <''(An,oA^i?)+ Yl W^{\n,2/k\R) 
n=l n=0 

= T2 + Ti + To" + . 

where the functions VFs(A//c^, i?) are obtained by collecting the contributions to (jl.ip coming 
from all the fields of spin s, and we have separated the h contribution from that of the other 
scalars (dubbed "np", for non-physical). Note that Ng is a function of R as well as of the spin s. 
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Explicitly, the functions Ws{z, R) in d = 4 are given by 



^ (z - 1) dtf'iR) - 2Rf"{R) + 2f'{R) - 2f'{R) 

W2 {z,R) = ^ . ^ ^ ^ , (5.5) 

{R-2)nR)-2f{R) 



W~^{z,R) 



Wi{z,R) = W^''{z,R) = -l, (5.6) 
?,{l-z^){wJ''{R)-mp){R)) 

(5.7) 



18f"{R)-2{R-3)f'iR) + Af{R) 

3(1 - z) (dJ'iR) - 2R}"{R) + 2/'(i?) 

+ 



lSf"{R)-2{R-2,)r{R)+^f{R) 

Note the simple form of Wi{z, R) and Wq^^^z, R), which is due to the ghost choice ()3.10p . With a 
standard ghost action we find slightly more complicated expressions, in particular with a depen- 
dence on 77 = —dt log Z, and for the fixed-point equation we have to require also r/ = —2, beside 
dtf{R) = 0. Once that is done, all the qualitative features of our analysis are left invariant. 

5.1 Approximating the sums 

Thanks to the cutoff choice (j5.2p . all the sums appearing in the FRGE can be performed analyti- 
cally. Indeed they all involve simple sums of the following type 

NsiR) — 



for m G {0,1,2}. 

Unfortunately, the use of (|5.2p is not safe from pitfalls. The main drawback is that since 
the upper bound on the summation, Ns{R), is a staircase- function of the curvature R, also 
the resulting sum Sm,s{R) is a staircase- function. Dealing with differential equations containing 
staircase-functions is quite unpleasant, and for this reason we will adopt some smoothing strategy. 
Of course such smoothing constitutes an approximation and it introduces an additional scheme 
dependence. We have explicitly verified that all the qualitative conclusions of our work are left 
unaltered by use of different smoothing choices. 

We now illustrate some possible smoothing choices, by means of the simplest example, that 
is, the function So^(){R) with uq = 0, also known as the spectral counting function. Specializing 
to d = 4 dimensions, we easily find 

So,o{R) = Vi[No{R)\) (5.9) 

where 

^ 1 

V{N) = Dn,o = -(1 + iV)(2 + Nfi3 + N) , (5.10) 



n=0 

and 



-m + \ RU8 + 25R) 
No{R) = ^-^ , (5.11) 
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Figure 1: A plot of the staircase function SQfi{R), together with four possible smoothing-curves, 
from top to bottom: the upper-edge, the asymptotic, the mean and the lower-edge interpolations. 

and where [.J denotes the floor function, which is what gives rise to the staircase nature of the 
function. 

In Ho] the simplest approximation was made, replacing [xj — )• x, that is, defining 

5iJ)(i?)=P(iVo(i?)). (5.12) 

This is the smoothing which we will call "upper-edge" , as it touches the original staircase function 
on all the upper edges of its steps. Analogously, we can define 

sj,'^\R)=V{NoiR)-l), (5.13) 

which we will call "lower-edge" , as it touches the lower edges of the steps (see Fig. [1]) . 

Such smoothing interpolations have the obvious disadvantage of containing square roots, which 
also are not particularly welcome in the differential equation. Rational interpolating functions 
can be obtained in several ways. One, is to consider the average of the upper- and lower-edge 
interpolations, which we will call "mean" interpolation: 

s!,~\r) = \ (5(J(i?) + 4:^{R)) . (5.14) 

We can also introduce an interpolation just by requiring that it matches the asymptotic leading 
behavior of the original function at i? — )• and at i? — t- -I-oo (a similar choice was used in |47] ) . 
We call such function the "asymptotic" interpolation, and in the specific example it reads 

5i7""'^^(i?) = ^ + l. (5.15) 

Finally, one could use the heat kernel interpolation, as tacitly done in j231l24j. On the relation 
between counting functions and heat kernel there is abundant literature (see for example [38] and 
references therein), and we will not indulge on that here. 

We found no evident qualitative difference between the mean, the asymptotic and the heat 
kernel interpolations, and we will report in the following only the results obtained within the 
asymptotic interpolation. 



15 



6 The fixed-point diff^erential equation 

The full FRGE in d = 4 dimensions takes the form 

^ [dJkiR) + mm - m'kiR)) =T2+Ti+ T^" + , (6.1) 

where we used (12. 5p and the explicit formula for the volume of S^, and where we have subdivided 
the r.h.s. into the contributions of the TT-tensor modes 



2^{dtnR)-2Rf"{R) + ^f'{R) 
T2 = — , (6.2) 



the vector modes 

the non-physical scalar modes (by which we mean all the scalars but the trace h) 

To' = -'-^^, (6.4) 

and finally the contribution of the trace mode h 

1 



2i?2 {-<df"{R) + {R- 3)f'iR) - 2f{R)) 

{ (^R" - 54i?2 _ 54) (dJ"{R) - 2Rf''\R)^ ^^'^^ 
- + I8i?2 + 12) (dtf'{R) - 2Rf"{R) + 2f'{R)) } . 

An analysis of such a partial differential equation seems a formidable task, and we will limit 
ourselves just to the associated fixed-point ordinary differential equation. 

At the fixed point dtf{R) = 0, and we can express the ODE in canonical form, by solving for 
/(3)(^), resulting m 

/(3)(^) = Af{fJ'J";R) ^ 

i?(i?4 - 54^2 _ f^^ _ 2)f'{R) - 2f{R)^ ' 



where 

AA(/, /', /"; R)=- 360Rf"{Rf + 768^^^ (^R^ -5R + 6^ f'{Rf 

+ f{Rf (^-276487r2/"(i?) + 30727r2(3i? - 5)f'{R) - 20R^ - 192^ 
+ f{R) (^-2 (^^^ + 18R^ + 45i?2 + 52R + 432^ f"{R) 
- 15367r2 (sR^ - WR + 6^ f'{Rf 

+f'{R) (276487r2(^ - l)f"{R) + 22R^ - UR^ + 192^ - 136j 
+ f'{Rf (-69127r2(^ - 2)Rf"{R) - 6^^ + 9^^ - 42^^ + 68i? + 216^ 
+ (r^ + 16^^ + 9R^ - 38R^ + 288R + 576) f'{R)f"{R) - QUATT^f{Rf 



(6.7) 
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6.1 Gaussian fixed point 

The Gaussian fixed-point solution is easily recovered as in [23]. We rescale f{R) ^ and 



look for a solution in the limit /3 — >• 0. The r.h.s. of (|6.ip is homogeneous of degree zero in / and 
its derivatives, while the l.h.s. is of degree one, hence in the limit of vanishing coupling we get 
the fixed-point equation 



interpretation is as usual: expanding the metric as Qf^i, = g^^ + \/]3H^^, where the background 
Qf^i, is a solution of the equations of motion, the only part surviving in the /3 — >• limit is the one 
quadratic in H^i,. 

In the following section we will investigate the possible existence of a non-Gaussian fixed-point 
solution of (|6.6p . 

7 Analysis of the fixed-point equation 

We are now going to discuss the following properties of the equation ()6.6p : fixed singularities at 
R = and R = movable singularities; asymptotic behavior at infinity. Performing numerical 
integrations and putting things together, we will draw some conclusions about the existence of 
non-trivial global solutions. 

7.1 Singularity at R = 0, and polynomial truncations 

As anticipated in Sec. El the fixed singularity at the origin is there for the same reason why the 
equation is of third order, i.e. because because f"{R) appears in the cutoff function ()4.8p and 
because of the dimensional nature of R. 

The presence of such a singularity acts in a way to reduce the number of independent initial 
conditions. This can be easily understood: a regular solution of an equation of the canonical 
type f"'{R) = F{f , f , f"; R) can be constructed in the neighborhood of any regular point Rq, 
substituting a series expansion f{R) = Z^„>oOn(^ — Ro)^, expanding in series the resulting F, 
and then imposing the equation order by order. The first of such equations will fix 03 as a function 
of oq, oi and 02, the next will fix 04, and so on, leaving the first three series coefficients free. 
However if F has a pole of order n at Rq, the r.h.s. will produce a Laurent series with n singular 
terms not matched on the l.h.s. (analytic by construction), and which have to be equated to zero; 
the first of such equations will only contain the lowest order expansions for /, /' and /", hence 
it will only depend on oq, ai and 02, thus providing a constraint on the initial conditions. This 
is precisely what happens at Rq = 0, where we have a simple pole in (|6.6p . 
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leading to the constraint (]2.9p . At higher orders, a^+s enters hnearly in Bn, hence the tower 
of equations obtained from (j7.ip can be solved iteratively, providing a series expansion of the 
solution, as function of the two initial data oq and Oi. 

The usefulness of the series solutions for the search of a global solution ends here, as in 
general it will only be valid at small R. In particular the convergence radius of the solution is 
generically limited by the presence of singularities in the complex plane, rendering it difficult 
to search for initial conditions leading to a regular solution on the whole real axis. We do not 
report here results obtained using either ratio or root test, as they performed very poorly. A 
more sophisticated analysis using Pade approximants was also attempted, but the presence of 
two initial conditions variables rather than one makes a systematic study a bit clumsy, and we 
gained from it no useful insights to report here. 

The polynomial truncations studied so far in the literature can be understood as an approx- 
imate way to reduce the oo'^ set of solutions to a finite subset by discarding the contribution of 
all the terms in the series expansion beyond a certain order A^. In other words, one makes the 
ansatz 

TV 

/(i?) = J^a^ii", (7.2) 

n=0 

and then imposes the first A^ + 1 equations coming from (|7.1|) . that is, two more than we would 
impose to solve at order A^ the Cauchy problem for (j6.6p with initial conditions /(O) = oq 
and /'(O) = oi. Yet another way to say it is that we solve ()7.ip up to aAr+2 (that is, up to 
Bn-i = 0), and then impose a^+i = aN+2 = 0. From the point of view of the differential 
equation, the imposition of such condition can be understood as an heuristic way to push farther 
away singularities [H]. 

Proceeding on such a route we find similar results to previous analyses. For example, trun- 
cating at A^ = 1 we recover the Einstein-Hilbert truncation, and we find a non-Gaussian fixed 
point at ao = 0.005984, ai = -0.01497, corresponding to A* = 0.1998 and G* = 1.329, with 
critical exponents 9± = 1.62425 it i3. 42642. Actually one finds a number of fixed-point solutions 
which increase with A^, but most of them are there for one truncation and not for another, or 
they are at unacceptable values (complex couplings, or negative Newton's constant). These are 
generally considered spurious fixed points and hence discarded, however a precise and reliable 
method to select which fixed points to keep and which to discard in this scheme is missing. In 
particular a systematic expansion and analysis of the convergence with A^ is very demanding from 
the computational point of view, essentially due the fact that the series solution of (j7.ip leads to 
high-order polynomials of two variables for the coefficients aN+i{ao-,ai) and aAr+2(ao, ai). 

7.2 Singularity at ^ = ^+ 

At -R = Rj^ = \J 3(9 -|- \/87) ~ 7.415 the equation has another fixed singularity in the form of a 
simple pole. The origin of such a singularity has to be looked for in the hh sector of the theory. 
With the operator choice in (|4.3p we have eliminated any singularity in the equation (j6.ip . but 
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we have shifted the eigenvalue of the constant h mode to a negative value: from Table [T] we have 
Ao,o = —R/3. Being negative, such mode is never excluded from the sum (the step function in the 
optimized cutoff only kills modes larger than k"^, which of course is positive). The f"'{R) term 
in (|5.7p is multiplied by a factor proportional to X]„>o(l — -^n. o/^^)' ^"^^ lowest mode 

this can become zero at some large value of R, where thus the equation develops a singularity. 
This is precisely what happens at 

As a result of the pole, if we Taylor expand the solution around i?+ we find again a condition 
that reduces the number of independent initial conditions to two. We write the expansion as 

f{R) = bo + bi{R- R+) + j{R- R+? + ^"(^ - ^+)" ' (^-3) 

n>3 

To leading order the equation (|6.6|) reduces to 

= ^^^ + o((i?-i?,)0), (7.4) 

R — R-\- ^ ' 

ij+(oo, 01,02) = , ^ = , (7.5) 

+ ^ ^ Vl^fmR\{hx{R+-2)-2hQ) ^ ' 



where 



and 



B+{hQMM) = - 768TT'^blR+ (rI - 5R++6) + bj ( -9^5|. + 366^^ - 68-R+ + 108 



+ 768TT^bl ({9b2 + 660)^+ - 2(962 + 1060)^+ + 126oj 

- 6162 (276487r26o(^+ - 1) + 63-Rij. + 826-R^ + 342i?+ + 144o) 

- 26061 (^15367r^6o(3^+ - 5) + URI - 7R\ + 96i?+ - 68 ] 
+ 26062 (^138247r26o + 18^^ + 99^^ + 52^+ + 486^] 
+ bl (614471^60 + 20^^ + 192) + 3606^-R+ . 



(7.6) 



In this case we find that the regularity condition is -B+(6o, 61, 62) = 0. 
7.3 Movable singularities 

It is well known that nonlinear ODEs can have also movable singularities beside the fixed ones. 

For our equation, we found that for generic initial conditions the solution ends at some R = Rc 
where the solution becomes singular. More precisely singular solutions can be constructed at any 
point R = Rc ^ 2, of the form 

fiR) ~ log \R - Rc\ (mo + mi{R - Rc) + O {{R - Rc)^)) + co + ci{R - Rc) + O ((i? - Rcf) , 

(7.7) 

where the only free parameters are co and Rc, all the others being determined as function of 
those. For example, 

Rl - 2-Rc - 54-Rj? + 108^2 _ 54^^ j_ 

mn = = • (7.8 

3456 7r2(fl,-2) ^ ' 
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7.4 Expansion at i? — oo 

The solution can also be worked out in the limit of i? — t- oo. This actuahy turns out to be the 
hmit in which the series expansion is most under control. 

We found that the asymptotic solution has an expansion of the following form: 



f{R)^Aie \ l + Y,dnR~ 

n>l 



(7.9) 



All the coefficients dn can be fully determined iteratively as function of only A, which is another 
crucial point for our results. The presence of only one free parameter in the asymptotic expansion 
is a very general property of the equation, deriving from the higher order nature of the pole at 
infinity (as perhaps more easily seen after mapping the point at infinity to a finite value, and 

We report here the expression for the first three terms 



studying the balance of terms in (j6.6 
of the asymptotic solution: 



d. 



-y, d. 



4(5 + 384^^2^ 
4 



= --(-17 + 17640A7r2 



1161216^^71^) 



(7.10) 



We have computed coefficients up to d^Q, and we found that in general the coefficient dn is a 
polynomial of order n—1 in A, with the coefficient of the highest power a number of order 10'^""^, 
suggesting that, at least for large enough A, the series should be convergent for R > i?< ~ lO^A. 
We have verified this conclusion by the root test (or Cauchy-Hadamard theorem), i.e. by a linear 
fi10 of Vn = \fd^ as function of 1/n to estimate ii< = lim„_).(x, r^. A plot of i2< as function of 
A is reported in Fig. [2j As expected, for sufficiently large values of A^ i?< grows linearly with 




-0.005 0.000 0.005 

Figure 2: A plot of the (inverse) radius of convergence i?< for the series in (|7.9p as a function of A (the 
series converges for R > i?<). The horizontal line is i? = 

^^We have checked that in the explored range of values for A the dependence of Vn on Xjn is quite close to linear, 
only with a slight convexity at certain small values of A. For comparison a ratio method was attempted too, i.e. 
taking r„ = |d„+i/d„|, but the resulting behavior seemed far from linear, and for this reason we stick here to the 
root test results. 
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A. At small A some spikes are observed in correspondence of the zeros of the coefficients dn{A). 
However, precisely at such values, r„ deviates sensibly from a linear dependence on 1/n, and the 
result from the root test are not reliable^ In any case we expect that in general the radius of 
convergence of the series will not extend beyond the fixed singularity at R^, with possibly the 
exception of accumulation point of the zeros of dn{A). 

An approach equivalent to the standard truncations can be applied also to the asymptotic 
series (|7.9|) . That is, a truncated fixed-point solution of order is obtained by solving the 
differential equation up to dN+i and subsequently imposing diy+i = 0. Of course we find in this 
way exactly N complex solutions for A. It turns out that there is a sequence ^i(A^) of zeros 
which converges very quickly to A* = —0.001663801, as shown in Fig. [3l while the others move 
in the complex plane in a characteristic pattern. Fig. HI 



-0.001663801461 - 
-0.001663801462 - 
-0.001663801463 - 
-0.001663801464 - 
-0.001663801465 - 
-0.001663801466 - 
-0.001663801467 - 
-0.001663801468 - 
-0.001663801469 - 
-0.001663801470 - 
-0.001663801471 - 



Figure 3: A plot of the most rapidly converging solution to dN+i = 0, as a function of iV S {1, 29}. 



Around a fixed point, we can evaluate the stability exponents 6i of the beta functions for the 
coupling of the asymptotic expansion!^ We have computed them for the fixed points of Fig. [3] 
up to order N = 8, and we report them in Table [2j It turns out that the number of relevant 
directions (corresponding to Re(0i) > 0) seems to grow with A^. 

We also computed the stability exponents for the other fixed points at yl < 0, which form 
another sequence A2{N) slowly drifting towards the first, and we report them in Table [3j Also in 
this case the number of relevant directions seems to grow with A^, with only one more irrelevant 



^^Note that at small values oi A the sequence r„ splits in two subsequences r2j and r2j+i, which differ significantly 
from each other at small j. However we can easily deal with such behavior by fitting appropriate subsequences 
of r„. On the contrary, in correspondence of the real part of the zeros of a coefficient dm{A), the corresponding 
element Vm will deviate from its neighbors, giving rise to an erratic behavior of the sequence r„(A) which is not 
easy to take care of. 

^■^Defining do = A, and dtdi — [3i{{d„}), the stability matrix is defined as Bij = djP^A* ■ The stability exponents 
{9i} are defined as minus the eigenvalues of Bij. 
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Figure 4: A table of plots of the structure of the fixed points for the large-i? truncations. The 
plots are ordered increasingly from top left to bottom right. For < 16, a positive real fixed 
point appears at odd values of the truncation order, and it is converted into couples of complex 
fixed points at even values of A. At A = 16 the real fixed point splits into a couple of real fixed 
point, one of which drifts towards the negative axis, and the other follows the previous patter, 
with odd and even orders inverted, until it undergoes a new splitting at A = 29. 
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N 


lO^Ai 


00 


01 


02 


03 


04 


05 


06 


07 


08 


1 


-1.319 


-4.26 


0.65 
















2 


-1.631 


29.2 


9.19 


-1.64 














3 


-1.660 


29.6 


7.47 


-1.80 


20.39 












4 


-1.663 


29.4 


7.25 


-1.82 


17.3-4.38i 


17.3-l-4.38i 










5 


-1.663 


29.7 


7.27 


-1.81 


14.6-4.24i 


14.6-|-4.24i 


20.9 








6 


-1.663 


29.5 


7.26 


-1.82 


14.0-3.19i 


14.0-1-3. 19i 


19.4-4.20i 


19.4-|-4.20i 






7 


-1.663 


29.8 


7.26 


-1.82 


14.4-2. 50i 


14.4-1-2. 50i 


16.7-4.72i 


16. 7-1-4. 72i 


22.5 




8 


-1.663 


29.5 


7.26 


-1.82 


14.8-3.55i 


14.8-1-3. 55i 


15.8-2. 32i 


15.8-l-2.32i 


21.6-l-3.83i 


21.6-3.83i 



Table 2: The critical exponents for the sequence ^i(A^) of Fig. [3l at the orders = 1, • • • ,8. 



N 


lO^Aa 


00 


01 


02 


03 


04 


05 


06 


07 


08 


3 


-0.225 


-19.1 


9.78 - 5.68i 


9.78 


-1- 5.68i 














4 


-0.522 


-27.3 


8.08 - 5.60i 


8.08 


-1- 5.60i 


-0.51 


13.0 










5 


-0.727 


-29.1 


7.68 - 5.96i 


7.68 


-1- 5.96i 


-0.51 


13.1 


14.9 








6 


-0.910 


-34.8 


8.55 - 6.05i 


8.55 


-1- 6.05i 


-0.56 


13.4 


15.0 -1- 1.75i 


15.0 - 1.75i 






7 


-1.05 


-47.8 


9.28 - 6.94i 


9.28 


-1- 6.94i 


-0.61 


14.2 


13.8 + 1.41i 


13.8 - 1.41i 


17.6 




8 


-1.17 


-78.1 


9.43 - 7.66i 


9.43 


-1- 7.66i 


-0.67 


13.5 


17.5 + 1.09i 


17.5 - 1.09i 


14.5 


15.3 



Table 3: The critical exponents for the sequence A2{N). 

direction compared to Ai{N), probably corresponding to the trajectory connecting the two fixed 
points. 

It should be mentioned that truncations with single terms added to the Einstein-Hilbert 
actions were considered in [23j . where it was observed that such new terms correspond to new 
relevant directions. Our results are hence compatible with that of [23]. 

Finally, it should also be noted that negative values of A correspond to an unbounded fixed- 
point action (see Sec. l7.6p . and might be not very appealing for this reason. Potentially interesting 
candidates for a fixed point could maybe be found also at positive A, slightly relaxing the type 
of convergence one is looking for. One way to visualize other interesting points is to show all the 
solutions for different values of N in the same plot, as in Fig. [5j In this way we can see that 
there are some accumulation points for the zeros of di\i-^-i{A), and these are potential candidates 
for global solutions of the differential equation. Of course if a sequence of zeros converges to a 
point A*, like in Fig. [3l this will be an accumulation point. On the other hand, accumulation 
points appear also whenever a fixed point which appears at odd orders, splits into a couple of 
nearby complex ones at even orders (or vice versa). One such sequence of points is present in our 
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3 


2.54 


-0.55 


6.60 


4.03 - 10.8i 


4.03 -1- 10.8i 










5 


4.65 


-0.52 


9.73 


2.12 - 18.4i 


2.12 -1- 18. 4i 


8.55 - 2.71i 


8.55 + 2.71i 






7 


5.72 


-0.54 


8.74 


9.54 - 4.15i 


9.54 -1- 4.15i 


12.7 


13.5 


-4.30 - 23. 7i 


-4.30 + 23. 7i 



Table 4: The critical exponents for the sequence of real fixed point at ^ > 0, appearing at odd 
values of N. 
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Figure 5: A simultaneous plot of all the fixed points for the large- i? truncations, at various N. 

solutions, and in Table U] we list the relative stability exponents. However, due to the fact that 
only odd orders appear, and that higher orders are very demanding from a computational point 
of view, the available data are not enough to establish any convergence of its critical exponents. 

7.5 Numerical integration 

One advantage of having a differential equation for the fixed-point effective action is the possibility 
of going beyond the series expansions discussed so far, and to study its solutions via numerical 
integration methods. In principle we would like to employ such methods in order to look as in 
|14| for the initial conditions leading to global solutions o Unfortunately, the presence of two 
free initial conditions at the origin, combined with complicated structure of the equation, and 
the presence of a fixed singularity at R±, renders such an analysis much more involved than in 
the scalar case. We report here on some preliminary results in this direction. 

The most efficient approach to the numerical investigations is to start integration at large 
R, where the series expansion presents only one free parameter!^ Imposing initial conditions 
matching (j7.9p . and integrating backward towards small R, we can monitor the appearance of 
singularities as we vary the parameter A. Naively we would like to vary A in order to find those 
values for which we can integrate all the way to i? = without running into a singularity of the 
type ()7.7p . However, because of the singularity at i?+ we expect that to be impossible, as indeed 
confirmed by explicit numerical integration. The plot Fig. [6] shows that the numerical integration 
breaks down once the fixed singularity at R^ is reached. As we already know from Sec. 17.21 in 

^*At least solutions defined on the whole positive axis plus the origin. As the equation (|6.6|l was derived assuming 
a spherical background, we might not trust it at negative R. 

^^The starting value of R has to be chosen within the convergence radius of the series in (|7.9|l . and also for this 
reason we have detailed the convergence analysis in the previous section. 
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Figure 6: A plot of the final integration point i?< reached numerically (starting at i? = 60), i.e. of the 
location of the movable singularity, as a function of A. The horizontal line is R = R+. 

order to have a regular solution at the condition i?_|_(6o, &i, ^2) = must be satisfied, which 
hence we would like to impose as a condition on A. However the numerical integration can only 
be carried out until an e away from We have thus evaluated and plotted B+{bo,bi,b2) at 
i?+ + e as a function of e for various values of A, in order to extrapolate the result to e = 0. It 
turns out that B+{bo,bi,b2) scales to zero with e in the range —0.0035 ^ A < 0.0005, while it 
diverges outside of that (if it reaches at all). As an example of scaling we report a plot of 
B+{bQ,bi,b2)/e as a function of e and A in Fig. [71 As a result, the regularity condition at R^ 



Figure 7: On the left we show a plot of i?+(&o,&i,&2)/e as a function of e at A = .00015 (similar 
plots are obtained in the whole range —0.0035 ^ A < 0.0005). As trivial as it might look like, the plot 
carries an important message: B+{bo, 61, 62) ^ e at these values of A. The plot on the right shows instead 
B-f-ibo, &i, 62)/e as a function of A: various plots for different values of e arc shown simultaneously, and the 
superposition is perfect. 

does not fix ^ to a discrete set but just to a continuous interval. 

Within such interval, numerical integration could be continued beyond ii^ by matching the 
numerical solution at R^ + e with the series solution ()7.3p . and use the latter to impose initial 
conditions at R+ — e for a new numerical integration towards R = 0. We would have to look for 
the range of A which allows integration down to R = e and here impose the analyticity constraint 
(j2.9p . again in the limit e — )■ 0. A number of technical difficulties appear along this path, in 
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particular due to a very slow convergence of the series ()7.3p . and a detailed implementation of 
this program is left for future work. 

However we have achieved an important result: the mismatch between number of conditions 
and number of parameters (two analyticity conditions for one free parameter) does not seem to 
provide an obstruction for the existence of global solutions. The condition at i? = leaves 
us still with a continuous degree of freedom, hence anything could happen when imposing the 
condition at i? = 0: we could find a discrete set (eventually empty) of solutions, as well as another 
continuous interval. 

7.6 Fixed point action 

One very general conclusion can be easily derived from the large R series expansion studied in 
Sec. 17.41 if ()6.6p admits a global fixed-point solution f*{R), then necessarily 

T* = TU^ = A* J d^x^R\ (7.11) 

for some finite A*. That is, the fixed-point effective action corresponds to an interacting R? 
theory. 

Here the effective action is the standard one, obtained from the average effective action in 
the limit /c — t- 0, i.e. it is the effective action obtained by a path integral with no infrared cutoff. 
Obtaining the effective action is generally a difficult task, which would require integrating the 
FRGE down to A; = 0, starting from some initial condition at some /cq > (e.g. at the UV cutoff, 
ideally with an initial condition on a trajectory emanating from the UV fixed point, to ensure 
independence from the UV cutoff). However, at the fixed point such integration is trivial, as the 
/c-dependence of T*^, is contained in a trivial scaling. We can easily prove (17. lip by noticing that 

n = k'' j d'^x^riR/k'), (7.12) 

and hence the limit /c — )• corresponds to the limit — )■ oo in f{R). Inserting (17. 9p . with A 
at the fixed-point value A* , into (j7.12p . and taking the limit A; — t- 0, we find ()7.1ip . Note that 
we had already met an R? action in the case of the Gaussian fixed point. However, in that case 
^* — 7- oo, requiring a rescaling of the fluctuation field, and a reduction to a free theory. 

The result (17. lip could be expected on dimensional grounds, from the simple fact that the 
fixed-point action should be scale invariant. However, naive scaling of a term of the Lagrangian 
is only valid if it has no anomalous dimension. For example, in d = 3 the scalar potential at 

6 

the Wilson-Fisher fixed point is V{4') = A*(j)^+'^ , and = in the LPA, but rj ~ 0.03 at higher 
orders of the derivative expansion |13j . What we found here is one more analogy between the 
f{R) approximation and the LPA: within the f{R) approximation the anomaouls dimension of 
R is zero and hence the fixed-point effective action is an R? action. 

Furthermore, such result was impossible to see in the polynomial truncations studied in pre- 
vious works: the results from such truncations \12\ [23l [24} [25] indicate that every coupling in the 
expansion is non-zero at the NGFP, and no hint on how the series would sum up was available. 
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8 Conclusions 

In this article we have derived and studied a differential equation for the renormalization group 
fixed points of gravity in the f(R) approximation. We have argued that such approximation plays 
a role analogous to the one played by the local potential approximation in scalar field theories, 
and taking seriously such perspective we have examined various properties of the solutions to our 
differential equation. In particular we have studied regularity conditions and large- i? expansion. 
Within the latter we have checked the convergence of fixed-point solutions obtained by truncations 
of the series expansion. 

The main results of our analysis are the following: 

• the fixed-point structure of large- i? truncations is much clearer than the one from small R 
truncations, and in particular we found striking convergence of a fixed-point solution (Fig. [3] 
and HI , which however turns out to have a growing number of relevant directions; 

• the analyticity condition at R = R^ is satisfied on a continuous interval of the parameter 
A, hence the presence of two fixed singularities in our equation is not an obstruction to the 
existence of global solutions; 

• if global solutions exist, they necessarily correspond to the effective action of an R^ theory. 

The main motivation for our investigations was the asymptotic safety scenario conjectured 
by Weinberg, and our hope is that the approach we have presented here will help understanding 
further the status of such conjecture. It should be pointed out however that our methods and 
results are quite generic and could also be useful for the investigation of other aspects of gravity, 
such as infrared modifications, stability of the de Sitter solution, and nucleation of black holes 
(see [l9] for a one- loop analysis of such questions within an f{R) setting). 

What we presented here was just a first step in the direction of exploiting the f{R) approx- 
imation in the spirit of the local potential approximation, and many more developments of this 
approach are possible and needed. For example, on a technical side, the use of a different cutoff 
would be desirable, possibly avoiding the staircase nature of the traces. Pushing the Table H] to 
higher orders would be important for testing the existence of a fixed point with A* > and with 
a finite number of relevant directions. And of course, a more comprehensive exploration of the 
numerical solutions is being performed, and it will be the subject of a future publication. 

Finally, the result of Sec. 17.61 suggests that in a more generic approximation (i.e. not relying 
on a maximally symmetric background) the fixed-point action would probably contain a Weyl- 
squared term. As a consequence, we feel that the issue of unitarity of the resulting theory will 
need to be discussed further, for example sharpening the argument presented in [26j . 
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9 Erratum 



Equation (j5.7p in the original paper contains a mistake, and it should be corrected to 

_ _ 3{l-z')(3dJ"iR)-6Rp){R)) 

W^{z,R) = , , ^ = ^ ^ 

18/"(i?)-2(i?-3)/'(i?) + 4/(i2) ^g^^ 

3(1 - z) [dJ'iR) - 2Rf"{R) + 2f'{R)) + 6 (/'(i?) + 6f"{R)) 

^ 18/"(i?) - 2{R - 3)f'{R) + if{R) ' 

The error propagates to Sections [6] and U\ affecting part of the numerical results, but not the 



qualitative aspects. In particular, equation (|6.5p should be 

I 

2i?2 (-9f"{R) + {R- 3)/'(i?) - 2f{R) 



R" - 54i?2 - 54] [dtf"{R) - 2Rp\R) j 2) 
- + 18i?2 + 12) (dJ'iR) - 2Rf"{R) + 2f'{R)) 
-36 {ie + 2) (/'(i?) + 6r(i?))}. 

Equations ()6.7p and (I7.6p should be corrected accordingly, as well as the numerical results in 
Sections 17.41 and I7.5[ All the qualitative aspects of the results (e.g. Section [7^ . as well as the 
general assessment and conclusions, are unaffected by such corrections. 
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